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I. INTRODUCTION 

In the Standard Model (SM) flavor-changing neutral 
currents (FCNC) are forbidden at tree level [1]. These 
decays are allowed at higher order and have been mea- 
sured in the K and B meson systems [2]. In the charm 
sector, however, the small mass difference between down- 
type quarks of the first two families translates to a large 
suppression at the loop level from the GIM mechanism 
[1]. To date, measurements of radiative decays of charm 
mesons are consistent with results of theoretical calcula- 
tions that include both short-distance and long-distance 
contributions and predict decay rates several orders of 
magnitude below the sensitivity of current experiments 
[3-8] . While these rates are small, it has been postulated 
that new physics (NP) processes can lead to significant 
enhancements [9]. 



In this paper we report results of a search for the FCNC 
decay of the neutral D meson into two photons. The 
only previous study was conducted by the CLEO col- 
laboration using 13.8 fb _1 [10]. Theoretical calculations 
predict that the decay D° — > 77 is dominated by long- 
distance effects. A calculation in the framework of Vector 
Meson Dominance (VMD) [11] yields 

B(^°^77) (VMD) -(3.5^)xl0- 8 . (1) 

A separate calculation using heavy quark effective theory 
combined with chiral perturbation theory (HQxPT) [12] 
reveals a similar dominance of long-distance over short- 
distance (SD) effects, with the SD branching ratio esti- 
mated to be [11] 

^ 77) (SD) * 3x 10- 11 . (2) 
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In the context of the Minimal Supersymmetric Standard 
Model, gluino exchange can enhance the SM rate by up 
to a factor of 200 [9] . The large number of charm decays 
in the BABAR dataset provide the opportunity to study 
this enhancement. A 200- fold increased rate would result 
in approximately 1370 events in the BABAR dataset (470.5 
fb" 1 ) with only six events predicted from the theoretical 
SM branching fraction (3.5 x 10 -8 , as determined in the 
VMD calculation.) A summary of the relevant branching 
fractions is shown in Table I. 

In this paper we also report a new measurement of the 
branching fraction for the decay D° —> tt°tt°, which is 
the dominant background in the D — > 77 analysis. 

TABLE I. Summary of predictions and measured values or 
limits for branching fractions relevant to this analysis. The 
results presented in this paper are not included in this table. 

Theoretical predictions 

Mode Value Reference 

D° -> 77 (SM,VMD) « (3.5 ±|"_g) x 1(T 8 Burdman [11] 

D° -» 77 (SM,HQ X PT) (1.0 ± 0.5) x 1(T 8 Fajfer [12] 

D° -> 77 (MSSM) 6 x 1(T 6 Prelovsek [9] 



Experimental results 



Mode 






Value 


Reference 


D° -> 


77 




< 2.9 x 10~ 5 


Coan [10] 


D° ->• 


7T 7T 


(2006) 


(7.9 ±0.8) x 10~ 4 


Rubin [13] 


D° — > 


7T 7T 


(2010) 


(8.1 ±0.5) x 10~ 4 


Mendez [14] 


D ^ 


KgTT 


(1.22 ±0.05) x 10" 2 


Nakamura [15] 



II. THE BABAR DETECTOR AND DATASET 

This analysis is based on a data sample corresponding 
to an integrated luminosity of 470.5 fb _1 collected by the 
BABAR detector at the SLAC PEP-II e+e~ asymmetric- 
energy collider operating at e + e~ center-of-mass (CM) 
energies of ^ = 10.58 GeV and 10.54 GeV. 

The BABAR detector is described in detail else- 
where [16]. Charged particle momenta and positions 
are measured with a five-layer double-sided silicon ver- 
tex tracker (SVT) and a 40-layer drift chamber (DCH). 
Charged hadron identification is provided by measure- 
ments of the ionization energy loss, dE/dx, in the track- 
ing system and the Cherenkov angle obtained from a ring- 
imaging Cherenkov detector (DIRC). An electromagnetic 
calorimeter (EMC) consisting of 6580 CsI(Tl) crystals 
measures the energy deposited by electrons and photons. 
These detector elements are located inside the cryostat of 
a superconducting solenoidal magnet, which provides a 
1.5 T magnetic field. The instrumented flux return (IFR) 
of the magnet allows discrimination of muons from pions. 

A detailed Monte Carlo simulation (MC) of the BABAR 
detector based on GEANT 4 [17] is used to validate the 



analysis and determine the reconstruction efficiencies. 
We use simulated events to optimize the selection criteria 
by maximizing significance, defined as N$ / y/N$ + Nb , 
where Ng and Ng denote the number of signal and 
background candidates in the MC simulation assuming a 
D° — > 77 branching fraction of 5.4 x 10~ 6 (five times less 
than the CLEO collaboration upper limit). The back- 
ground samples include e + e~ — > cc, e + e~ — > qq, q — u, d 
or s, e + e~ — > B Q B . and e+e - — > B + B~ decay modes. 

III. EVENT RECONSTRUCTION AND 
SELECTION 

For the decay modes used in this study we require 
that the neutral D meson originates in the decay D* + —> 
D°tt + , which is referred to as a D* + tag. (The inclusion 
of the charge conjugate modes is implied unless otherwise 
stated.) Without such a tag the signal is dominated by 
combinatoric background. To avoid uncertainties in the 
number of D* + mesons in the BABAR dataset, we per- 
form measurements of the D° — > 77 and D° — > tt°tt q 
branching fractions relative to a well measured reference 
mode. The D° — > K s ir° decay is chosen for this purpose 
due to its large branching fraction of (1.22 ± 0.05)% [15] 
and partial cancellation of systematic uncertainties in our 
measurements. 

In the D° — > 77 analysis the D° candidate is formed 
by combining pairs of photon candidates. D° candidates 
are required to have an invariant mass between 1.7 and 
2.1 GeV/ c 2 for the D° — s- 77 analysis and an invariant 
mass between 1.65 and 2.05GeV/c 2 for the D° —> tt°tt° 
analysis. A photon candidate is defined as energy de- 
posited in the EMC, which is not associated with the 
trajectory of any charged track and which exhibits the 
appropriate shower characteristics with a lateral moment 
[18] greater than 0.001. The photon candidates are se- 
lected to have CM energies between 0.74 and 4 GeV. 

In the D° — > ir Q ir° analysis two tt° candidates each with 
CM momentum above 0.6 GeV/c are combined to form a 
D° candidate. The ir° candidates are formed by combin- 
ing two photon candidates with lateral moment less than 
0.8. The list of ir° candidates also includes single EMC 
clusters containing two adjacent photons (merged tt°). 

The D° candidates for the D° — > K^tt reference 
mode are formed by combining a tt° candidate as de- 
fined above with a Kg candidate consistent with the de- 
cay Kg — > 7r + 7r~. The 7r + 7r~ invariant mass is required 
to be between 0.491 and 0.505 GeV/c 2 . To be selected 
as a K s candidate the decay length significance must 
be greater than 3, where the decay length significance 
is defined as the measured flight length divided by its 
estimated uncertainty. 

In all modes, D° candidates are combined with tt + 
candidates selected from tracks with CM momentum be- 
tween 0.05 and 0.45GeV/c. A kinematic fit is applied to 
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the events, requiring the candidate D° invariant mass to 
be between 1.6 and 2.1 GeV/c 2 . Both the D° and 7r + arc 
constrained to originate from a common vertex within 
the beamspot to satisfy the D* + tag requirement. 



IV. BACKGROUND STUDIES 

Backgrounds from B meson decays are removed by 
selecting D* + candidates with CM momentum greater 
than 2.85 GeV/c in the case of D —> 77 and greater than 
2.4GeV/c in the case of D° -> 7r°7r°. The difference re- 
flects cuts optimized to separate MC samples. 

In order to minimize systematic uncertainties the ref- 
erence mode analysis was performed separately for each 
of the two signal modes, each time using identical crite- 
ria that were optimized for the respective signal mode. 
These selections result in 95% rejection of B meson de- 
cay modes. The D° —¥ 77 decay mode has significant 
backgrounds due to QED processes, which are largely 
removed by requiring that the total number of charged 
tracks in the event be greater than four and the number 
of neutral candidates in the event be greater than four. 

The dominant background to D° — > 77 is due to D° — > 
7r°7r° decays. To remove this background, we implement 
a 7T° veto. From our sample of D° — > 77 candidates 
we reject all events in which one of the photons can be 
combined with any other photon candidate in the event 
to form a ir°. This veto rejects 95% of the background 
and keeps 66% of the signal. 

The D° — > 77 analysis signal efficiency is 6.1% with 
the corresponding reference mode (D° — > KgTr°,Kg — > 
tt + it~) efficiency at 7.6%. The D° — > ir°ir° analysis sig- 
nal and reference mode efficiencies are 15.2% and 12.0%, 
respectively. 



V. FIT PROCEDURE AND RESULTS 

For each of the three decay modes we determine the sig- 
nal yield using unbinned maximum likelihood fits to the 
invariant mass distribution of D° candidates passing the 
above selection criteria. The overall probability distribu- 
tion functions (PDFs) are sums of functions describing 
signal and background distributions obtained from the 
Monte Carlo simulation. The relative normalizations of 
these functions are free parameters while the individual 
shapes are fixed. 

In the D° — > 77 analysis the signal PDF consists of 
a Crystal Ball [19] function and a bifurcated Gaussian 
distribution. The background PDF is a 2 nd -order Cheby- 
chev polynomial and the D° — > 7r°7r° background shape 
is described by a second Crystal Ball function. In the 
D° —> 7r°7r° analysis the signal is described by a sum 
of a Gaussian, a bifurcated Gaussian, and a Crystal Ball 



function, and a background PDF described by a 3 rd -order 
Chcbychcv polynomial. 

The invariant 77 mass distribution obtained from the 
£)° — > ryry analysis is shown in Fig. 1 together with pro- 
jections of the likelihood fit and the individual signal and 
background combinations. The signal yield is — 6 ± 15, 
consistent with no D° —> 77 events. We convert this re- 
sult to a branching fraction for D° — > 77 relative to the 
D° — > Kgir reference mode using 



B(D° -> 77) = 



-N(D° -> 77) 



-N(D° -> K° s tt°) 



xB(D° -> K° s tt°), 
(3) 

where N and e are the yield and efficiency of the re- 
spective modes and B(D { 



is the known D° — > 



t + tt branching fraction [15]. In this anal- 
ysis the D° -> K%n° signal yield is 126599 ± 568 events. 
We find B{D° -> 77) = (-0.49 ± 1.23 ± 0.02) x 10" 6 
where the errors are the statistical uncertainty and the 
uncertainty in the reference mode branching fraction, re- 
spectively. 
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FIG. 1. The 77 mass distribution for D° — > 77 candidates 
in data (data points). The curves show the result of an un- 
binned maximum likelihood fit to the measured mass distri- 
bution. The solid blue curve corresponds to signal compo- 
nent resulting in a slight negative yield, the long-dash red 
curve corresponds to combinatoric background component, 
and the small-dash pink curve corresponds to the combina- 
toric background plus D° — ► 7r°7r° background shape. The 
X 2 value is determined from binned data and is provided as 
a goodness-of-fit measure. The pull distribution shows differ- 
ences between the data and the solid blue curve with values 
and errors normalized. 



The invariant mass distribution for events in the D° — > 
7r°7r° analysis is shown in Fig. 2. The signal yield for 
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D° -> ir°ir° is 26010 ± 304 events. For D° -> K%tt° 
(mass distribution not shown) the signal yield is 207538 
± 1143 events. Adjusting Eq. 3 for the D° — > 7r°7r° case 
we convert this yield to a branching fraction and find 
B{D° -> 7r°7r°) = (8.4±0.1±0.3) x 10~ 4 . The first error 
denotes the statistical uncertainty and the second error 
reflects the uncertainties in the reference mode branching 
fraction. 
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dates in data (data points). The curves show the result of the 
unbinned maximum likelihood fit to the measured mass dis- 
tribution. The solid blue curve corresponds to the full PDF 
including the signal and the dashed red curve corresponds to 
the combinatoric background component. The \ 2 value is de- 
termined from binned data and is provided as a goodness-of-fit 
measure. The pull distribution shows differences between the 
data and the solid blue curve with values and errors normal- 
ized. 



VI. SYSTEMATIC UNCERTAINTIES 

Several systematic uncertainties cancel partially or 
completely when the branching fraction is measured with 
respect to the D° —> Kgir reference mode. The uncer- 
tainty in tracking efficiency and vertcxing 1.39%. The 
uncertainty due to photon reconstruction efficiency in the 
ratio of the signal mode branching fraction to the refer- 
ence mode branching fraction is 3.0% and 0.6% for the 
D° — > tt°tt° and D° — > 77 analyses, respectively. 

In order to account for the uncertainty arising from 
fixed PDF shapes, the parameters determined from the 
Monte Carlo simulation, are varied by random amounts 
sampled from the covariance matrix retaining correla- 
tions among parameters. The values of these parameters 
are fixed and the resulting PDF is fit to data allowing the 



yield to float, the lcr width of the obtained signal yield 
distribution is taken as the systematic uncertainty. In the 
D° — > 7r°7r° analysis, fixing the signal and combinatoric 
background shapes results in 0.20% and 0.80% systematic 
uncertainties, respectively. Fixing the D° — > K^tt signal 
and background shapes for the reference mode results in 
0.17% and 0.63% systematic uncertainties, respectively. 

Potential differences in ir° veto efficiencies between 
data and the Monte Carlo simulation are estimated using 
a sample of candidates for the physically forbidden decay 
D Kgj. The difference in the ratios of numbers of 
candidates before and after the veto between data and 
the Monte Carlo simulation is taken as the systematic 
uncertainty. We measure the difference as a function of 
the number of photons in the event and as a function of 
the photon energy. In all cases, the variations are found 
to be less than or equal to 1.8%. 

In order to account for imperfect modeling of D* + 
hadronization, a 4% correction is applied to the MC for 
normalized momenta, x = p{D* + )/p ma _ x (D* + ), within 
the region x = 0.575 to x = 0.7 to match cross- 
section measurements made by the CLEO collabora- 
tion [20]. Wc calculate the ratios of signal efficien- 
cies (ffDO-). 77 /ff£)0_ ) .K0 7r o ! £_D°->7r 7r°/ e -D -s.if°7r°) wrt h an d 

without this correction applied to the MC and deter- 
mine systematic uncertainties of 0.02% and 0.03% for 
the D° —J- 77 and D° —> 7r°7r° modes, respectively, due 
to this correction. 

To account for systematic uncertainties due to apply- 
ing a particular set of selection criteria, we vary the selec- 
tion criteria and recalculate the results. The reconstruc- 
tion efficiency is determined from MC and the efficiency- 
corrected yield is measured from data when each set of 
selection criteria is applied. These yields are found to be 
distributed normally and the standard deviation is taken 
to be the systematic uncertainty. Choosing particular 
event selections for the D° — > 7r 7r° and D° — > Kgir 
studies results in systematic uncertainties of 2.50% and 
0.76%, respectively. 

The systematic uncertainties are summarized in Table 
II. For the D° — > 7r°7r° mode a total systematic uncer- 
tainty of 4.2% is obtained by adding all contributions in 
quadrature. 

For the D° —> 77 analysis we combine all system- 
atic uncertainties with the statistical uncertainties in 
the upper-limit calculation. In a Monte Carlo simula- 
tion study we generate event samples using the complete 
background PDF from the data fit and repeat the branch- 
ing fraction calculation 14000 times varying all sources of 
systematic uncertainties in the process. For each branch- 
ing fraction calculation the selection values on the contin- 
uous variables are varied within ranges established from 
the D° — s- 7r°7r° analysis. In each calculation the pa- 
rameters of the signal and background PDFs are varied 
within their uncertainties while fully accounting for the 
correlations among them. Systematic uncertainties such 
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TABLE II. Summary of systematic uncertainties, a, for each 
mode. [*] These sources of systematic uncertainty are assessed 
in a combined Monte Carlo simulation study. See text for 
details. 



Source of a(D° ->■ 77) a(D° -t ir°w°) 



Systematic Uncertanity (%) (%) 

Tracking (K° s ) and Vertexing 1.39 1.39 

Photon Reconstruction 0.60 3.0 

tt° Veto 1.8 

D* + Hadronization 0.02 0.03 

Signal Shape * 0.20 

Background Shape * 0.80 

Selection Criteria * 2.5 

D° -> A'gvr Signal Shape 0.10 0.17 

D° -» K^ir Background Shape 0.53 0.63 

D° -> K° s ir° Selection Criteria 0.76 0.76 

Total Systematic Uncertainty * 4.3 



as tracking, photon reconstruction, it veto, and D* + 
hadronization are added in quadrature and the D° — > 77 
signal efficiency is varied randomly according to a normal 
distribution with a width equal to this total uncertainty. 
To account for the uncertainty in the branching fraction 
of the D° — > KgTT reference mode, the nominal value is 
varied randomly according to a normal distribution with 
a width equal to the established D° — > Kgir branching 
fraction uncertainty [15]. The resulting distribution of 
B(D° — > 77) branching fractions is shown in Fig. 3. 



ID 




B(D°^yy) 

FIG. 3. Distribution of branching fraction calculations when 
varying all sources of uncertainty. 

Integrating this distribution to 90% of the area with 
B(D° — > 77) > gives us the expected sensitivity of our 
analysis. We find 

B(D° -> 77) < 2.4 x 1CT 6 (4) 



at 90% confidence level. 

VII. RESULTS 

In this paper we present a new measurement for the 
branching fraction for a D° meson decaying to two neu- 
tral pions: 

B(D° -> 7r°7r°) = (8.4 ± 0.1 ± 0.4 ± 0.3) x 10~ 4 , (5) 

where the errors denote the statistical, systematic, and 
reference mode branching fraction uncertainties, respec- 
tively. 

We also report the result of a search for the decay of a 
neutral D meson to two photons. The observed yield is 
— 6± 15 consistent with no D° — > 77 events. Our analysis 
has an expected sensitivity of B(D a —> 77) < 2.4 x 10 -6 
at 90% confidence level. In order to obtain an upper 
limit for B(D° — > 77) we repeat the sensitivity study 
described in Section VI with the signal yield set to the 
measured value of —6 instead of and find 

B(D° -> 77) < 2.2 x 10" 6 (6) 

at 90% confidence level. 

This result is consistent with our expected sensitivity 
and with SM expectations. As stated earlier, gluino ex- 
change has been postulated to possibly enhance the SM 
rate by up to a factor of 200 [9]. Based on the upper limit 
of the branching fraction for D° — > 77 presented in this 
paper, the enhancement of the rate over the expected SM 
rate cannot exceed a factor of 70. 
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